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-5 ■ ABSTRACT 

6 

' We study the cosmological evolution of scalar fields that arise from a phase transition at some energy 

C/3 i scale Ac. We focus on negative power potentials given by = cA^+"(^~" and restrict the cosmological 

^ ' viable values of Ac and n. We make a complete analysis of V and impose SNla conditions on the dif- 

ferent cosmological parameters. The cosmological observations ruled out models where the scalar field 
, has reached its attractor solution. For models where this is not the case, the analytic approximated 

5-H ' solutions are not good enough to determine whether a specific model is phenomenologically viable or 

not and the full differential equations must be numerically solved. The results are not fine tuned since 
a change of 45% on the initial conditions does not spoil the final results. We also determine the values 
of iVc, Nf that give a condensation scale Ac consistent with gauge coupling unification, leaving only 
four models that satisfy unification and SNla constraints. 
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1 INTRODUCTION 



Recent cosmological results given by the superonovae project SNla and the Maxima and Boomerang 
observations have lead to conclude that the universe is flat and it is expanding with an accelerating 
velocity. This conclusions show that the universe is now dominated by a energy density with negative 
pressure with fip/j, = 0.7 ± .1 and < —2/3 Q . This energy is generically called the cosmological 
constant. An interesting parametrization of this energy density is in terms of a scalar field with 
gravitationally interaction only called quintessence Q . The evolution of scalar field has been widely 
studied and some general approach con be found in jll]. The evolution of the scalar field 
depends on the functional form of its potential V and a late time accelerating universe constraints the 
form of the potential pl| . 

In this paper we will concentrate on negative power potentials because they lead to an acceptable 
phenomenology and because they are naturally obtained from gauge group dynamics. Negative power 
potential have been extensively studied first by and then as tracker fields by [Q. Steinhardt 

et al. 1^ showed that scalar field with a negative power potential V = c(f>~"' with n > 5 has already 
reached its tracker solution but are not cosmologically acceptable because they have > —0.52. 
However, if the scalar field has not reached its tracker solution by today, we will show that the models 
may lead to an acceptable phenomenology and the final results depend on the initial conditions and 
on the value of n. Contrary to the tracker models, no analytic solution is good enough to determine 
the value of w^o (from now on the subscript o will refer to present days values) and it is sensitive to 
the whole dynamics. We will solve numerically the differential equations and we will constrain the 
values n of cosmologically viable models, including the big bang nucleosynthesis (NS) constraints [Q. 
We will also give approximated analytically solutions. 

Tracker solutions are widely favored because they do not have a fine tuning problem of initial con- 
ditions. But even more, they are independent of the initial conditions since the range of the initial 
conditions can vary up to 100 orders of magnitude. The models with n < 5 do depend on the initial 
conditions but it is important to remark that they do not have a fine tuning problem. The initial 
conditions can vary up to 45%, the solutions are still fine and the values of the initial conditions are 
completely "natural" , i.e. they are of the same order of magnitude as the other relevant cosmological 
parameters. So, to conclude, one thing is to have a model with no dependence on the initial conditions 
and another is to have a fine tuning problem. "Natural" models in physics should not have fine tuning 
problems but they do in general depend on the initial values as it is the case for our models. However, 
for any initial conditions we will end up with —1 < < — 2/(2 + n) = wtr, where n gives the inverse 
power and wtr is the tracker value. 

Negative power potentials Q]-^ can be obtained using Affleck et al. ("ADS") superpotential [p^ . 
The condensation scale Ac of the gauge group SU{Nc) can be determined from a high energy scale (A) 
using renormalization group equation in terms of Nc, Nf and it is then naturally to ask if it is possible 
to have a common gauge coupling unification with the Standard Model ("S.M.") gauge groups 
We will show that this is indeed possible and we will give the values of Nc,Nf where gauge coupling 
unification is achieved. 

The cosmological picture in the case of gauge coupling unification is very pleasing. We would assume 
gauge coupling unification at a scale A for all gauge groups (as predicted by string theory) and then 
let all fields evolve. At the beginning all fields, "S.M." and those from the SU{Nc) gauge group, are 
massless and redshift as radiation until we reach the condensation scale Ac. Below this scale the fields 
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of SU{Nc) group will dynamically condense and we use ADS potential to study its cosmologically 
evolution. Interesting enough, the relative energy density of the SU{Nc) group fl^ drops quickly, 
independently of the initial conditions, and it is close to zero for a long period of time, which may 
include nucleosynthesis, until very recently (around one e-fold of inflation). The energy density of 
the universe is at present time dominated by the scalar field with ~ 0.7 and a negative pressure 
< —2/3 leading to an accelerating universe 

The paper is organized as follows. In section |^ we give the general frame work to derive the scalar 
potential for (j) from non-abelian gauge dynamics using ADS potential. In section |^ we analyze the 
cosmological evolution of (p and we concentrate in the non attractor regime. We derive analytic 
formulas for and (jjo as a function of the initial conditions and on n and we discuss in detail the 
possible choices of initial conditions and show that the models do not have a fine tuning problem. In 
section ^ we constrain the values of n in order to have w^o < —2/3 while in section ^ we comment 
on the possibility to have models with a gauge coupling constant unified with the couplings of the 
standard model and we explicitly give these models. In section ^ we give further examples. Finally, 
we summarize and conclude in section m. 



2 POTENTIALS OF THE FORM V = cAt 



In this work we study quintessence field (scalar) cf) with negative power potentials that arise from a 
phase transition at some stage of the evolution of the universe. The energy scale of the phase transition 
is given by Ac and the initial value of the scalar field (j) is naturally given by (pi = Ac since it is the 
relevant scale for the transition. The potential we will consider is of the type 

V = c2A^+'V~" (1) 

with c a constant (we will comment on the value of c in section ^ and has a runaway behaviour. 

This class of models would have a vanishing potential V = for energy scales above Ac since the 
phase transition has not taken place yet and there is no (p field. At the phase transition energy scale 
Ac a potential l^(Ac) ~ A^ and a field 0(Ac) = Ac are generated. Below Ac, the (j) field becomes 
dynamically and it evolves to its minimum. The cosmological evolution depends on the functional 
form of V and for eq. (|l]) with n > we expect (p to roll down its potential. This class of potential has 
been chosen because they can be obtained from a phase transition of non-abelian gauge dynamics (see 



subsection 2.1) ||15| and because they lead to a quintessence interpretation of the (p field. However, if 
other physical process leads also to inverse power scalar potential the cosmological evolution studied 
in section ^ and the conclusions remain valid. 

The energy scale Ac is expected to be considerably smaller than the reduced Planck Mass, mpi, 
so the initial value = is much smaller than 1 and this has interesting consequences for the 
cosmological evolution of (we will set from now on the reduced Planck mass to one mp; — G/Stt = 1). 

The normalization of the field is important and we will consider, for simplicity, the p field to be 
canonically normalized, Lk = (K^)~^\dfj.(p\'^ with — 1 and = dK/dcp. However, the complete 
Kahler potential K is in general not known. The canonically normalized field <p' can be defined by 
<p' — g{(p, (j))(p with the function g given from solving the differential equation iiT^ = (5 + + 4'g,j/f' ■ 
For ^ 1 we do not expect any large contributions to the kinetic term but for 4> 1 the Kahler 
potential could give a significant contribution and could spoil the runaway and quintessence behaviour 
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of (f). In order to see this, we can expand the Kahler potential as a series power K ^ |0P + Siai|(/)p/2i 
with Oi some constants of order one and to be determined by the specific model. If we approximate, 
for simplicity, the canonically normalized field 4>' by 4>' = The potential in eq.(|l]) would 

then be given by 

= c2Ar"r"(A'^)t-i (2) 

For n <2 the exponent term of K'jj^ in eq. is negative so it will not spoil the runaway behaviour of 
<j) but for 71 > 2 the extra terms could stabilize the potential. In the lack of a better understanding 
of K we will work with canonically normalized fields but we should keep in mind that for n < 2 the 
results are robust while for n > 2 the contribution from the Kahler potential could spoil our results 
and must be determined. 

2.1 ADS Potential 

The potential in eq. (|l|) can be obtained from the non-perturbative dynamics of a non-abelian asymp- 
totically free gauge group SU{Nc) with Nf chiral + antichiral fields Q in N=l supersymmetric theory. 
At energy scales much larger than the condensation scale the gauge coupling constant is small and 
the Q fields are free elementary fields. As the universe expands and cools down, the energy of the 
elementary fields Q becomes smaller while the gauge coupling constant grows. When the gauge cou- 
pling constant has the critical value to condense the Q fields, then all the elementary fields will no 
longer be free and they will form "mesons" and "baryons", as in QCD. This effect takes place at the 
condensation scale Ac and below this scale the correct description of the dynamics of the non-abelian 
gauge group is in terms of the condensates 4>. In order to study the dynamics of these fields 4> we 
use the ADS superpotential, which is exact (i.e. it does not receive radiative or non-perturbative 

contributions) and it is given by W{(l)) — {Nc — Nf) ( j^^^QQ^ j 0- In terms of the gauge 

singlet combination of chiral and antichiral bilinear term (f> =< QQ > the globally supersymmetric 
scalar potential is given by eq.(||) with c — 2Nf, det < QQ >= n^j(/)| and n = 2 + f/^'p^^ 0"|9)- 

If we wish to study models with < n < 2, which are cosmologically favored as we will see in section 
we need to consider the possibility that not all Nf condensates cj)i become dynamically but only a 
fraction v are (with Nf > > 1) and we also need Nf > Nc- It is important to point out that even 
though it has been argued that for Nf > Nc there is no non-perturbative superpotential W generated 
Ip^ this is not always the case . The possibility of having ^ Nf can be done with gauge group 
with unmatching number of chiral and anti-chiral fields or if some of the chiral fields are also charged 
under another gauge group. In this later case we have c = 2^, n = 2 + '^-N/ ^'^^ ~ ^ condensates 
fixed at their v.e.v. < QQ >= @- 

3 Evolution of 

We will now determine the cosmological evolution of a scalar field </> with an inverse power potential 
regardless from what physical process it arose. We will concentrate on potentials that have no reached 
the tracker solution yet since they can give the correct value of w^o and we will give approximately 
analytic solutions to w^o and (/>o- 
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The cosmological evolution of with an arbitrary potential V{(j)) can be determined from a system of 
differential equations describing a spatially flat Friedmann-Robertson- Walker universe in the presence 
of a barotropic fluid energy density that can be either radiation or matter, are 

H = -^{p-,+Pj + ^^), 

p = -'SH{p + p), (3) 

where H is the Hubble parameter, </> = d(j)/dt, p {p) is the total energy density (pressure). We use the 
change of variables x = and y = and equations (||) take the following form |l^ : 

2-^2/^ + -^x[2x'^ + 7^(f - x'^ - y^)] 
VN = -^Xxy+^y[2x^ + j-,il-x^ -y^)] (4) 
Hn = -^H[2x^ +j^il-x^ -y^)] 



1 



where N is the logarithm of the scale factor a, N = ln{a); /at = df /dN for / = cc, y, iJ; 7^ = 1 

and A(-/V) = —V'/V with V' = dV/d(j)- In terms oi x,y the energy density parameter is fi^ = a;^ + y^ 

2 _ 2 

while the equation of state parameter is given by W(Ij = p^/p^ = frijip-- 

The Friedmann or constraint equation for a flat universe fl-y + fi^ = 1 must supplement equations (|^) 
which are valid for any scalar potential as long as the interaction between the scalar field and matter 
or radiation is gravitational only. This set of differential equations is non-linear and for most cases has 
no analytical solutions. A general analysis for arbitrary potentials is performed in the conclusion 
there is that all model dependence falls on two quantities: X{N) and the constant parameter 7^. In 
the particular case given by y cx 1/0" we find A in the asymptotic limit. If we think the scalar 
field appears well after Planck times we have Xi = n mpi /Ac 3> 1 (the subscript i corresponds to the 
initial value of a quantity). An interesting general property of these models is the presence of a many 
e-folds scaling period in which A is practically a constant and Sl^ ^ 1. Figure |l| shows the rapid arrival 
and long permanence of this parameter to its constant value, together with the final decay to zero. On 
this last regime we have A ^ that implies — < and ^ > leaving us with Qj; = x^ ^-y"^ ^ 1 

X y 1 — V 

2 _ 2 

and Wrfio = ^2_^y2 ^ — 1, which are in accordance with a universe dominated by a quintessence field 
whose equation of state parameter agrees with positively accelerated expansion. The development of 
fl^ can be in agreement with the restriction of the nucleosynthesis stage 570(NS) < 0.1 iQ as well as 
with the observational result Q^o = 0.7 (the subscript o refers to present day quantities). This can be 
observed in Figure |^, together with the evolution of Wff, which fulfills the condition < —2/3 

The analysis of inverse power potentials has been extensively studied 0-||ll|- However, the analysis 
has not been specific enough to determine their viability to describe the late evolution of our universe. 
In Steinhardt et al. Q the scalar field was required to track before present day and this imposes a 
constraint on n to be larger than 5 and thus ruling this models out since they have w^o > —0.52 in 
contradiction to SNla data. The models we will concentrate on are, therefore, models with n < 5 
where A has not reached its tracker value. 
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Figure 1: Evolution of A for n = 1. The vertical line marks the time at fl^o — 0.7 with Ntotai — 17.04. 
For future reference we give now the scaling value ofn (j) M 



(t>sc = (t>i + V^^^i ^<pi < 

= 0. + x/6 (^-i= + ^^05[y-^]) for r!^, > 1/2 (5) 

The scaling value only depends on the initial conditions ll^i, it is independent on Hi, since 0^ ^ 1. 
The tracker value of w is given by 

fi 

Wtr = -I + J^iw^ + 1) (6) 

and it is an attractor solution valid for large n, when cj) is already tracking. In the tracker limit Q, 
i.e. n — 5, from eq.(H) one has wtr = —0.28 but the value obtained numerically is only Wo — —0.52 
for > 0.25. For smaller n the discrepancy is even worse since the scalar field has not reached its 
tracker value, obtained from eqs.(0) and ([l5|), 



'?i(2 + n)il0o 

"tr = Y ^ ' (') 

which is larger than (j)sc if f^^j < "^^j^^"^ ^^o. 

A semi-analytic approach is useful to study some properties of the differential equation system given 
by eqs.(|4[). To do this we initially consider only the terms that are proportional to A, since A 1, then 
we follow the evolution of x, y and H so every period has a characteristic set of simplified differential 
equations. The parameter fi^ is adequate to divide the process into four periods, the first one being 
a short lapse in which — constant, easier to recognize in Figure ^; the second is defined from the 
fall of this parameter to negligible values; the third is the so-called scaling period and in the fourth 
is again considerable, eventually reaching the value 0.7. 

The phase plane x — y provides an illustrative approach, useful for the analysis. We see, from Figure 
^, that the system follows, at first, a circular path with = x"^ + y"^ = fl^i constant and ends up with 
~ fl^i ^ y^. Then, x and y decrease to negligible values, this situation prevails all along the 

Our value of (l>ac differs in the case of n^,, > 1/2 by a factor of 1/^2 from § and they use Mpi = 1 
instead of nip = Mpi/y/Sn = 1 as we do 
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Figure 2: Evolution of fl^ (solid curve) and (dashed curve) for n = 1. The vertical line marks 
the time at fi^o — 0.7. 

scaling period. Finally, a growth in both parameters causes ^Ic/, to reach what we set as the final value 
0.7, preferred by observational results. From the restriction over the equation of state parameter, 
< —2/3 and the observational range for Q^o = 0.7 ± 0.1, we can define a region limited by the 
expressions: — jj^^x^ with = —2/3 and = fi^o — with 0.6 < il^o < 0.8. 
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Figure 3: Phase plane x — y for n = \. Starting point is (0,0.5). The region defined by < —2/3 
and il^o — 0.7 ± 0.1 is shown. 

The minimal value j/mm of y after its initial steep descent is given from eq.(^ with yjv = 0, ~ il^^ 
and \xy/ij2 = -Hn/H = 37^/2 by 

A, 



(t>min — — -\/6r20i (8) 

and we have approximated Hmm — Hi in eq. (^ . Shortly after y reaches its minimum value the scaling 
period begins. In this period we neglect the term proportional to A in eqs.(^ to find: 

y ' N 
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which leads to yH = Hminymin — Hiymin- Notice that the dependence of A from y and H is given 
hy X = A (j/iJ)^/" with A a constant, therefore from eq.(|^) we have A = constant (i.e. (f) = constant) 
during aU of the scahng period, this holds for any n. Furthermore, we may neglect squared terms on 
X and y in the third equation of system (||), since they are small, to get the expressions 

y = ymme*-^-^ (10) 

The quantity y has an increasing exponencial form for almost all of the process, so the duration of 
this regime can be seen as the total time (see Figure Now, in order to calculate the number of 
e-folds from the initial value to present day we consider equation (|l^) to end up with: 

Ntotal = lH^] (11) 

with ymin given by eq.(||) and yo — 0.8 (to have Vl^o — -7,11)^0 < —2/3. The evolution of ln(.T), ln(y) 
and hi{H) as a function of A'' is seen on Figure (^). 
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Figure 4: Evolution for n = 1 oiln{x), ln{y) and ln{ilcf,) (dashed, dotted and solid lines, respectively). 
The vertical line marks the time at VIaio = 0.7. 



If we consider eqs.(lO) and assume that the end of the scaling period is very close to today we get 
an approximated equation yoHo ~ Hminymin = Hiyi{4>i/ ipmin)- This, together with eq.(^ and the 
definition of y^ = c^A'^/3Hf establishes an expression for Ac, the energy scale at which the scalar field 
appears, in terms of yo, Hq-. 

Ac=f^^V"</'^ (12) 



and (j)min = f'vSn^i. The later expression is a semi-analytic calculation of the initial energy scale 
of a specific model. Finally, the value yo is set to be of the order of 0.8 to satisfy simultaneously 
^if>o — 0.7 and condition w^o < —2/3. 

Of course, we could have guessed expression (^ using the definition of j/^ — — 5^2— 0o " to give 

A^.f^^Mii^^.H* (13) 
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and the last equality holds approximately since we expect to lia,ve yo , (j^o 

of the order of one. 

Now, we wish to determine the value of w^o, Vo and (/>o. We use the differential equation for = + 1 
and IITtI 



(j^)^ = 37^(2- 70) (^A^^-lj 

{n4,)N = 3(7^-7^)17^(1-17^). (14) 

We see that 7^ is extremized at it;^ = 7^ — 1 = —1, 1 and at = —I + ^ = — 1 + ^3^. We have 
checked that the value of at the maximum evaluated at (jjo is a very good approximation (within 
5% of the numerical value) 



+ (15) 



Eq. (15) should be compared withW which differs by the factor of Qr^o — 0.7. Of course, if we do 
not have the exact value of (j>o eq.([l5|) is not very useful. In order to determine (/)o we evolve eqs.(p^ 
and eq. (jj) from present day values to the scaling regime where w^sc — — 1 and x'^^ <^ yl^ ~ ^4>sc with 
the condition Vt^o ~ 0.7. This evolution is model independent if < + 1 <C 1 |l^ and from the 
definition of we have 

^A4+" = ylHlc^- ^ ylHl^-, (16) 



and we obtain 

2 tt2 

VoVo - -0sc^'0o (17) 



where we have used i/Ar/i/ ~ -3{l-n^)/2, H^J = {l^n^o)e^^^ and -17^ = ^,{,06-^^^ /{I - 
+ r!^oe-3AJV)^ AiV » 1 0. Using eqs.{^, and (|) we can solve easily for 0o and/or yo 
in terms of IJ^o, n and Q^i (via 0sc), 



or equivalently 



6(O0o - yl) 



y 1^00 = (18) 



i:,170o. (19) 



In order to analytically solve eqs.dlq,tt9) we need to fix the value of n and we can determine w^o by 
putting the solution of ( [T^ ) into eq.(|l5|). Eq.(|l8|) can be rewritten as (po — 0sc(l — "^f70o/6(/)o)^^/" and 
we see that > 4>8c and 0o > «\A7^io/6 and we see that (j)o is of the order of 1 (fi^o ~ -7) regardless 
of the initial conditions. However, the exact value does indeed depend on the initial conditions but 
for any initial conditions we will have — 1 < w^o < wtr- 

If l<t>o — 'n?^4,ol^<t>t ^ 1 '^^^'^ Ii^s — 4>sc and for the simple cases of 7i = 1,2 and 4 we can solve 
explicitly for and we find 0o|„=i = 0sc/2+^9(?!)2^ + 6^4,0/6, yl\n=i = (/'sc(-3(/'sc+ V^'/'L + Sil^o, ) 

(po\n=2 = v/^irT2fV73, yl\,,=2 = 3</)Lf^0o/(30L+2f^0o) and(/)o|„=4 = y^417^o/3 + ^Qf^f, + l&n^o/'^, 
yl\n=4 = ^4,0 — 8^7^^/(417^0 + Y^90f^ + IGO^o), respectively. Notice that the value of (l)o,w^o at 
17^0 = 0.7 does not depend on Hi or Ho and it only depends on fl^i (through ^^c) and n. 

We show in Figure]^ how varies for different initial conditions 17^^ with n = 18/7 ~ 2.57 fixed. 
We see that for larger fJ^j we end up with a smaller lu^o and this is a generic result as can be seen 



8 



from eqs.(|l^),(p[) and ( |15| ) since for larger fi^^ one has a larger (/)sc, (t>o and therefore a smaller w^o- 
This can be seen also from Figure |^ where for a small Q^i a plateau arises in (the field (f) has 
already reached its tracker value by present day). From eqs.(p^),(p^ we notice that for smaller n one 
gets a smaller w^o (see Figure |^) . 



3.1 Initial Conditions 

It is well known that inverse power law potential lead to attractor solutions and are therefore 
independent of the initial conditions (one has 100 orders of magnitud range for Q^i). However, as 
pointed out by Steinhardt et al, this is only the case for n > 5, for smaller values of n the field (p 
has not necessarily reached its tracker value at present time. However, even though in these cases the 
present day quantities depend on the initial conditions there is no fine tuning problem since we can 
vary the initial conditions on a wide range of values (i.e. 45%) and the end results are still physical 
acceptable. 

The differential equations given in eqs.(^ depend on values of x,y,X = n/(j) but they do not depend 
on the absolute value of H, i.e. we have the same evolution for x, y, H as for x, y, H' = kH where 
k is an arbitrary constant. This scaling freedom allow us to set the normalization of H as we wish 
and in particular to have il0o = 0.7 at = Ho for any values of the initial conditions Xi,yi. This 
implies that we can have a quintessence model for arbitrary initial conditions. Once Ho is fixed we 
get Ac ~ i/o^*-^^"-* (eq.(pO|)) and Hi — cA^/j/i (from the definition of yi). One of the main difference 
with the tracker solution is that or (j)o do not have the same values in all cases but depend on fJ^^ 
as can be seen from eqs.([l5|) and (|8|). However, the dependence on il^^ is mild since for example for 
n = 1 one has that (j)o only depends as y/Tl^. For this reason once Ho and are fixed we still have 
a wide range of values of Q^i giving the correct phenomenology, i.e. Afl^i = 45%. For <C 1, i.e. 
yi = cAc/Hi ^ 1, with a fixed Ac one has a larger Hi and the time of expansion up to present day 
is also larger. In this case it is posible that the (j) field has already reached its tracker value, even for 



n < 5, and we would have w^o — Wtr = —2/{2 + n). Indeed, it can be seen from eq.(18) that a Q^i <C 1 
renders & (po n-^/fi^o/S and independent of the initial conditions^. For any initial conditions we will 
end up with — 1 < w^o < Wtr- 



4 Quintessence Restriction on n 

Before analyzing the quintessence restriction imposed on n we would like to comment on the value of 
c in the potential eq.(|l]). It can take different values with different physical interpretations. 

For simplicity of arguments let us assume that we have no kinetic energy at the beginning and that 
p<t>i = Vi = c^Ac. The initial energy density is fi^i = p^ij Pc- Taking p^i = c^Ac, with c = and 
Pc = we have that the initial energy density is given by fi^i — {2v)'^ [Kc / . In this case we see 
that the initial energy density depends on the ratio between the condensation scale Ac and the critical 
energy density E = pV^. 

Another possibility is to take p^i = g^iT"^ and pc — grotT'^ where grot, g4> are the total and quintessence 
number of degrees of freedom, respectively, at a temperature T. In this case we have Vt^i — g^/grot 

*We would like to stress out that eqs.(|l8|) and (^^ are valid only for < + 1 <C 1 and in the tracker 
regime this region is no longer satisfied. 
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Figure 5: Variations on fi^i lead to different physical situations given by w^o- We have taken 
ri^i = 0.9,0.25, 10"^" dashed, solid and dotted lines, respectively. The vertical line marks the time at 

n^o = 0.7. 
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Figure 6: Evolution of Vl^ and for n — 1. Parameter 17^ evolves in a similar way for fl^i = 0.25 
(solid line) and fi^i = 10^^" (long-dashed line), varying only in the total time, Ntotai = 16.7 or 21.1 
in the previous order. A discrepancy between the two cases at fi^o = -7 is seen for parameter 
because it has, for = 0.25 a local maximum with w^o = —0.93 (dotted line) while for fi^^ = 10"^" 
(short-dashed line) a plateau appears with iii^o = —0.76. The vertical lines mark the time at fJ^o = 0.7. 
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and it only depends on the ratio of quintessence and total degrees of freedom and not on the energy 
scales. The condensation scale is then Ac = {g^j (?YI^T and p^i ~ g,pT^ = c^A* = Vi. 

Let us now, study the restrictions to the values of n. In both cases, mentioned above, we get similar 
restrictions so we will only consider the first one. 

For a fixed n we take the following conditions: ho must be in the range .7± .1 and Q^o must belong to 
the interval 0.6 — 0.8. We can restrict n according to different physical arguments. These restrictions 
are encountered while solving numerically the differential system settled by equations (^), the first 
restriction comes from the observational value of the relevant parameter w^o- The limit w^o < —2/3 
can be translated to w^ff < —0.7 Notice that in this potentials w^jf < w^o contrary to the 
general arguments and the reason is that is still growing by today. From this analysis we find 
that, in order to fulfill the Weff condition, n must be smaller than 2.74 as shown in Figure for 
= 0.25. The value of w^o depends not only on n but also on fi^i and it decreases with increasing 
fJ^i. If we fix n and increase fi^i we find that examples originally discarded by the Weff < —0.7 
restriction now enter the physically permitted group. The example with n = 3 is depicted in Figure 
^. An equipartition value of Qr^i is 0.25 but if we allow il^i to be as large as 0.75 than the restriction 
on n is only n < 5.2. 

-0. 65 

-0.7 
-0.75 
Weff -0.8 
-0.85 

-0.9 
-0. 95 

1 1.5 2 2.5 3 

n 




Figure 7: Restriction over n from Weff < —0.7 with c = 1 and fi^i = 0.35) 

Other restriction comes from big bang nucleosynthesis (NS) results that require f2^(NS) < 0.1 at the 
energy scale range of NS: 0.1 — 10 MeV [Q. To account for this we have to either consider il^j < 0.1 
or, for example, if fJ^i = 0.25 we must take out 1.2 < n < 2.1 because for this range the initial value 
Hi lies within the range of values of H^g. For all the values of n allowed by Weff and NS restrictions, 
a variation of 45% on Q^i can be performed without disturbing the permanence into the observed 
ranges of Hg and Q^o- 



5 Unification of Gauge Couplings Constants 

The condensation scale Ac used in eq.(|l|) is from an elementary particle point of view an arbitrary 
scale that sets the energy scale of the phase transition. However, if the inverse power potential eq.(||) 
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-0.725 
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-0.775 
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Figure 8: Restriction We// < —0.7 avoided by increasing with n — 3 {c— 1) 

is obtained from a non-abelian asymptotically free gauge group then we can relate the condensation 
scale Ac to other energy scales using the renormalization group equation. The one-loop evolution of 
the gauge coupling constant for an SU{Nc) gauge group with Nf chiral fields gives a condensation 
scale (strong coupling g^^ ^ 1) 

Ac = A,e"^ (20) 

where bo = (3iVc - iV/)/167r2 is the one- loop beta function and A^^gx are arbitrary energy scale and 
coupling constant, respectively, which include high energy scale where the original chiral fields Q are 
weakly coupled. 

It is well known that the gauge coupling constants of the standard model get unified at a energy scale 
Agut — lO^^GeV with a coupling constant ggut — \/ 47r/25.7. 

We want to impose to our model gauge coupling unification, i.e. the coupling of the gauge group 
responsible for quintessence should be unified at Agut with the standard model gauge groups j9|. In 
this case we require A^; , g^ in eq. (|2^) to take the values A^ — Agut , dx — Qgut and we have A^ — 
Agutexp[—1 /2boggut]- This is not a necessary condition but opens the possibility of thinking of the 
model as coming from string theory after compactifying the extra dimensions or on a Grand Unification 
Scheme where all gauge coupling constants are unified. 

Of course, not all values of Nc, Nf will give an acceptable phenomenology. This is because the 
cosmological evolution of (p and the gauge coupling unification set independent constrains on the 
condensation scale Ac and on Nc,Nf. From a cosmological point of view Ac depends on the inverse 
power n (see eq.(|l^)), which is a function of Nc,Nf, and from gauge unification Ac depends also on 
Nc, Nf through bo (see eq. (^) ) . These two constrains reduce drastically the allowed values of Nc, Nf 
and we also require Nf > 1 and iVc > 2 to be integers. 

In table Q we give the different values of A^c, Nf, v for which we have gauge coupling unification. We 
can see that there are only a small number of possible models (11). The first five models have an N f 
which differs form an integer by less than 0.05 while the other 6 models differ at most by 0.10. All 
other combinations of A^c, N f, v have a larger discrepancy and do not lead to Ac = A„. If we further 
constrain the models to agree with the cosmological observations (i.e. w^o < —0.7 requiring n < 5) 
we are left with only 4 models (number 1,2,3,11 of table |l|). All of these 4 models have n <2 and the 
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Num 


Nc 




V 


n 


Ac(GeV) 


1 


3 


5.98 


1 


0.66 


6 X 10"^ 


2 


6 


14.97 


3 


0.66 


6.9 X 10"*^ 


3 


7 


18.05 


4 


0.55 


1.6 X 10-s 


4 


8 


5.97 


5.97 


13.83 


1.3 X 10^2 


5 


8 


6.96 


3 


13.55 


7.6 X 10^^ 


6 


3 


1.90 


1 


5.66 


1.3 X 10^ 


7 


5 


3.91 


2 


9.38 


4.7 X 10^ 


8 


6 


5.09 


2 


10.85 


3.7 X 10^0 


9 


7 


5.08 


5.08 


12.64 


3.9 X 10^1 


10 


8 


20.90 


4 


6.75 


2.4 X lO"'^ 


11 


8 


21.10 


5 


0.47 


5.7 X 10"^ 



Table 1: Models with = A„. The first 5 have A^/ < 0.05 while the other 6 have AA^/ < 0.10 
discrepancy from an integer 

quantum corrections to the Kahler potential are, therefore, not dangerous. Notice as well that only 
two models (4,9) have v — Nf and in both cases n > 12. 

As an example of a model with gauge coupling unification we have a gauge group Nc = 3 with A^^ = 6 
and u ~ 1 which entails a value n = 2/3 according to n = 2 + ^ ^-Nf ~ ^/"^ 0- -'^'-"^ ^^^^ model 
we find, from the numerical solution, a total time Ntotai = 12.96 which does not superposes on the 
NS range 19.6 — 27.2. The Weff < —0.7 restriction is satisfied with w^ff — —0.98 and the condition 
from experimental central values ho — 0.7 and fi^o = 0.7 is also fulfilled taking fi^j = 0.25. The 
condensation scale is Ac = 4.2 x 10~* GeV. A full analysis of this model is presented in ]T^ . 

6 Further Examples 

We have already given some example in the previous sections. In table |^ and ^ we show the numerical 
results for different values of n with initial condition Cl^i ~ 0.25 and for different initial condition with 
n=3 fixed, respectively. 

Other interesting examples are when — Nf. The condition n < 2.74 (for fl^i < 0.25) requires 
Nf/Nc < 0.15 and therefore Nc > 7. For Nc = 8,Nf = 1 one has n = 18/7 ~ 2.57 and using 
fl^i = 0.25 and one obtains w^o — — 0.65, w^ff = — O.73,0o = 2.01, A^o = 30.3 and A/, = 2.3 GeV 
Ac = 1.6 GeV while for Nc^7,Nf = 1 one has n = 8/3, w^o = -0.64, Wef / = -0.71, (po = 2.04, No = 
30.9, Mi = 4:.2GeV, Ac — 3GeV. However, these models do not have Ac — A„, i.e. they are not unified 
with the SM gauge groups. 
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7 Conclusions 



We studied negative power potentials and we constrain the initial conditions an the power of the 
potential to satisfy the SNla results. For n larger than 5, the scalar field cf) has already reached its 
tracker value and is too large. So, we need to concentrate on potentials with n < 5 to comply 
with SNla results. We gave a semi analytic solution to w^o and in terms of Ho,n^i and n and 
we have solved numerically for some relevant cases. We obtained that w^o depends on Sl^^jn and it 
decreases with increasing fl^i while it becomes smaller for larger n. If we assume equipartition initial 
conditions with il^i < 0.25 than n is constrained to be smaller than n < 2.74, however, if we allow 
for fl^i — 0.75 the constraint is relaxed to n < 5.2. We have shown that one can vary the initial 
conditions up to 45% without spoiling the observational cosmological values at present time. For any 
initial conditions we will end up with — 1 < w^o < Wtr- 

We have seen that the negative power potentials can be derived from Affleck's potential and in order 
to avoid problems with the Kahler potential one requires n < 2 which implies that Nf > Nc and 
that not all condensate become dynamically (i.e. ly ^ Nf). For v — Nf one needs Nf/Nc < 0.15 to 
have w^o < —2/3. Furthermore, we have shown that it is possible to have a quintessence model with 
gauge coupling unification for all gauge groups, Standard Model and the gauge group responsible for 
quintessence, but the number of models is quite limited (4). 

This work was supported in part by CONACYT project 32415-E and DGAPA, UNAM project IN- 
110200. 
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n 


Hi (GeV) 




K (GcV) 




Weff 


Ntotal 


1/2 


1.33 X 10" 


ss 


7.48 X 10^9 


3.74 X 10' 


-9 


-0.97 


-0.98 


10.72 


2/3 


1.16 X 10" 


-33 


6.99 X 10"^ 


3.49 X 10" 


-8 


-0.97 


-0.98 


12.96 


1 


3.46 X 10- 


-30 


3.82 X 10-6 


1.91 X 10" 


-6 


-0.93 


-0.95 


16.97 


2 


4.68 X 10- 


-22 


4.44 X 10-2 


2.22 X 10- 


-2 


-0.77 


-0.83 


26.33 


18/7 


1.68 X IQ- 


-18 


2.66 


1.33 




-0.65 


-0.73 


30.42 


3 


3.41 X lo- 


-16 


37.92 


18.77 




-0.57 


-0.65 


33.07 



Table 2: Numerical solutions for different values of n with c = 1 and f2^, = 0.25. 





H, (GeV) 


1/4 


Ae (GeV) 




Weff 


^total 


1 X 10-^ 


5.48 X 10- 


-14 


480.90 


19.12 


-0.55 


-0.50 


35.69 


0.25 


3.34 X 10- 


-16 


37.54 


18.77 


-0.57 


-0.65 


33.07 


0.5 


2.85 X 10- 


-16 


34.68 


20.62 


-0.76 


-0.83 


32.88 


0.75 


3.05 X 10" 


-16 


35.88 


23.61 


-0.88 


-0.92 


32.75 


0.99 


5.13 X 10- 


-16 


46.53 


32.82 


-0.97 


-0.98 


32.20 



Table 3: Numerical solutions for different values of fi^j with c = 1 and n = 3. 
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